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ABSTRACT 


Luckring, James Michael. A Method for Computing the Core Flow in 
Three-Dimensional Leading-Edge Vortices (under the direction of 
Dr. Fred R. DeJarnette.) 

A theory is presented for calculating the flow in the core of a 
separation-induced leading-edge vortex. The method is based on matching 
inner and outer representations of the vortex. The inner model of the 
vortex is based on the quasicyl indrical Navier-Stokes equations; the 
flow is assumed to be steady, axially symmetric, and incompressible and, 
in addition, gradients in the radial direction are assumed to be much 
larger than gradients in the axial direction. The outer model is based 
on the three-dimensional free-vortex-sheet theory, a higher-order panel 
method which solves the Prandtl-Glauert equation including nonlinear 
boundary conditions pertinent to the concentrated vorticity 
representation of the leading edge vortex. The initial conditions for 
the inner representation draw upon conical flow notions and include an 
asymptotic viscous subcore. All matching is based on identifying 
overlapping zones among the various theories which are either derived 
analytically from asymptotic concepts or demonstrated numerically. The 
resultant flow is evaluated a posteriori for evidence of incipient 
vortex breakdown and the critical helix angle concept, in conjunction 
with an adverse longitudinal pressure gradient, is found to correlate 
well with the occurrence of vortex breakdown at the trailing edge of 
delta, arrow, and diamond wings. 
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NOMENCLATURE 


All nonreference quantities are nondimensional, unless otherwise 
indicated. Lengths are referenced to wing root chord, velocities to 
the free-stream reference velocity, and pressures to twice the free- 
stream dynamic pressure. Nondimensional ization is explicitly shown 
for deviations from this convention or for purposes of clarity. 
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1.0 INTRODUCTION 


Interest in vortical flows embraces a rather broad spectrum of 
aerodynamics and, for that matter, fluid mechanics. Aerodynamically, 
vortices are a natural consequence of lifting flow and are manifested 
at scales ranging from a boundary layer thickness to characteristic 
wing dimensions. They can persist for miles downstream from a typical 
transport aircraft or, on the other hand, undergo an abrupt disruption 
commonly referred to as bursting or vortex breakdown. In addition, 
vorticies are a usual consequence of flow instability, play a role in 
the large scale structure of turbulence, and are a dominant feature of 
separated flows. They are prevalent in flows such as a draining basin 
and can be manifested on geophysical scales in the form of tornadoes 
and hurricanes. The fluid mechanics of vortical flows has been 
addressed by several authors, among them Kuchemann (1966) and, in a 
broader sense, lugt (1983). 

The research of one class of vortical flow, the separation 
induced leading edge vortex phenomena, is largely motivated by the 
advent of high speed aircraft with supersonic capability. The design 
of these aircraft for efficient supersonic cruise often requires 
consideration of slender wing concepts which include comparatively 
high leading-edge sweep angles, thin wing sections, and in some cases, 
relatively sharp leading edges. At the higher angles of attack 
encountered for take-off and landing as well as for subsonic-transonic 
maneuver, the subject vortex flow naturally occurs on these wings and 
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can be exploited for a variety of beneficial effects. The application 
of slender wing benefits to military aircraft has recently been 
reviewed by Polhamus (1983). 

The computation of these vortex flow effects which are dominant 
for aerodynamic applications has proven to be a challenging task. 

Even for simple three-dimensional geometries with sharp edges at 
subsonic speeds, it has only been during the last ten years or so that 
methods have been developed to provide reasonable estimates of 
inviscid force-moment properties and, subsequently, pressure 
distributions. Most of these methods are based upon explicit 
representations of the vortex by various approximate means, as opposed 
to solving a globally applicable formulation such as the Reynolds 
averaged Navier Stokes equations or, possibly, the Euler equations. 
Apart from the usual drawback of turbulence closure for the Navier- 
Stokes equations, problems affiliated with large memory requirements 
along with the large execution times for these equations are 
particularly exacerbated for the subject flow. As opposed to attached 
flows where large gradients occur in narrow regions near the wing . 
surface, the zone of high gradients for the vortex flow case can 
easily extend for distances on the order of the local wing semi span or 
more, thus requiring a comparatively large zone over which a fine grid 
would have to be employed. This situation is particularly acute when 
viewed in the context of current limitations to algorithms for solving 
three-dimensional, nonlinear flows as discussed by South (1985). 

While progress in this area will no doubt be realized, the presently 
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available methods with explicit vortex modeling will continue to be 
useful for some time. This will particularly be true if they are 
extended by embedding representations of additional flow effects (such 
as those due to viscosity), much as has been done for attached flows. 

The purpose of the present research is to formulate a composite 
representation of a fully three-dimensional, separation-induced, 
leading-edge vortex with emphasis on the flow in the vortex core. The 
primary objectives of the research are, first, to calculate the 
general properties of the vortex core flow for the fully three- 
dimensional outer flow of simple wing shapes, and, second, to examine 
the resultant flow a posteriori for evidence of incipient vortex 
breakdown. The approach taken is to embed an inner representation of 
the core flow in an outer representation of the overall wing/vortex 
flow field. With this approach, considerable advantage can be taken 
of the extensive modeling studies that have been performed for 
isolated vortex cores as well as for three-dimensional leading-edge 
vortices, so long as appropriate matching conditions can be 
established between the two models. 

The present investigation draws heavily on the isolated vortex 
core solutions of Stewartson and Hall (1963) and Hall (1967) as well 
as the vortex sheet solution of Johnson, et al. (1980). These 
theories were found to be appropriate for the present study and 
offered adequate growth potential for future investigations. The 
composite representation of the vortex is formulated by matching an 
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inner representation of the vortex core which accounts for viscosity 
and continuously distributed vorticity to an outer representation of 
the vortex which is inviscid and has concentrated sheets of spiral 
vorticity. A sketch illustrating the theoretical model of the present 
investigation is presented in figure 1.1. 

Select background and review material is presented in section 
2. Included are a discussion of prominent leading-edge vortex flow 
features as well as a focused review of vortex core theory, leading- 
edge vortex theory, and vortex breakdown theory. Some additional 
review material addressing these disciplines in a broader sense is 
also cited. The particular theoretical methods, as implemented in the 
present formulation, are then developed in detail. In section 3, the 
initial plane solution for the inner representation of the vortex core 
is presented. The solution is conical and includes a viscous sub core 
which is asymptotic to an analytical solution of the incompressible 
and conical Euler equations. This solution is advanced in space with 
the nonconical quasicylindrical Navier-Stokes equations, section 4. 

The three-dimensional boundary conditions for the inner vortex core 
representation are provided by the free-vortex-sheet theory (section 
5) and the matching of these solutions is described in section 6. 
Results and discussion of the coupled formulation are presented in 
section 7. Future considerations based upon these results are 
contained in section 8, and some general conclusions are presented in 
section 9. Some frequently used parameters are shown in figure 1.2. 
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2.0 BACKGROUND AND REVIEW 

2.1 Leading-Edge Vortex Flow Features 

The basic flow structure of the separation-induced leading-edge 
vortex is shown in figure 2.1. For slender wings with relatively 
sharp leading edges, the flow is dominated at moderate to high angles 
of attack by spiral shear layers of concentrated vorticity emanating 
from the leading edges. These shear layers roll up in the presence of 
the wing under their own influence into what is commonly referred to 
as the vortex core; collectively, this flow structure constitutes the 
primary vortex. The primary vortex induces a transverse boundary- 
layer flow on the wing upper surface which can separate and cause a 
smaller secondary vortex. The chief attribute of the overall flow is 
the development of large vortex lift increments which are nonlinear 
and stable at moderate to high angles of attack. The stable vortices 
persist well downstream of the generating wing and eventually diffuse. 

As angle of attack is increased further, the stable vortex flow 
pattern drastically changes as the vortices exhibit an abrupt 
breakdown (or burst) downstream of the generating surface. With 
additional angle-of-attack increases, the burst point propagates 
rapidly upstream, and over the lifting surface. A water-tunnel 
photograph from Erickson (1981) illustrates strake vortices bursting 
in the vicinity of the wing trailing edge and empennage for a YF-17 
model (figure 2.2). Although the present effort is directed primarily 
at delta wings, this photograph emphasizes the complexity of the 
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vortex breakdown flow field for practical geometries. The occurrence 
of breakdown in the vicinity of the lifting surface has many adverse 
effects; these include lift loss, pitchup, buffet onset, control 
decrements, nonlinearities, and configurational sensitivities. As a 
consequence, the "usable lift" range for configurations employing the 
vortex lift concept is often limited to angles of attack sufficiently 
below the critical angle of attack (ab) at which bursting has reached 
the near field of lifting or control surfaces. The near-field 
breakdown stability boundary, ab, is presented in figure 2.3, along 
with the near-field vortex asymmetry boundary, aa» Tor sharp edged 
delta wings at incompressible conditions. This figure (after Polhamus 
(1971)) is a summary of experimental data; there are presently no 
theoretical methods available for predicting these boundaries or 
aerodynamic properties under conditions when the angle of attack 
exceeds the appropriate critical value. 

Some examples of the adverse effects affiliated with near-field 
vortex breakdown are presented in figures 2.4 and 2.5 for a 65° delta 
wing. In figure 2.4 the water-tunnel photograph of Lambourne and 
Bryer (1962) illustrates clearly the flow of a core streamline near 
the axis of a vortex undergoing breakdown; both axi symmetric and 
spiral type core flows are demonstrated directly down stream of the 
burst point. Also shown, the lift losses due to breakdown are 
substantial. The theoretical lift from which experimental losses are 
judged was determined with the free-vortex-sheet theory of Johnson, et 
al. (1980). For this case, the losses occur gradually and precede 



7 


ab. At lower angles of attack, some of the differences between the 
theory and experiment are also due to the relative accuracies of both. 

The effects of breakdown on longitudinal stability and buffet 
properties are equally pronounced as shown in figure 2.5. The 
experimental pitching moment data of Wentz and Kohlman (1968) 
demonstrate a marked pitchup somewhat prior to the onset of 
breakdown. The correlation between theory and experiment at low to 
moderate values of lift coefficient is less accurate than usually 
achieved with the FVS theory, primarily because of experimental 
effects such as free-stream curvature and support-strut inter- 
ference. (This aspect of the pitching moment data has been confirmed 
in a private communication with Dr. Wentz.) The dynamic data of 
Boyden and Johnson (1982) exhibit buffet onset in the vicinity of ab 
as indicated by the dynamic root bending moment coefficient. In this 
experiment, the buffet intensity for a > ab was chiefly a function of 
reduced frequency, although buffet onset itself was insensitive to 
this parameter as well as to Reynolds number, dynamic pressure, and 
subsonic Mach number. 

Other aerodynamic consequences of near-field breakdown can be 
equally pronounced, particularly with regard to lateral-directional 
effects as discussed by J. L. Johnson, et al (1980). Collectively, 
these adverse effects emphasize the need to be able to predict this 
phenomena. 
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2.2 Vortex Core Theory 

■ 1 i 

Representations of the flow in the core of a vortex generally 
include the effects of distributed vorticity and viscosity. Virtually 
all of the theories are axi symmetric, and most model the vortex in 
isolation from natural surroundings for simplified boundary 
conditions, e.g., the models are usually appropriate for a vortex in a 
tube, or at best, conical edge flow. Newman (1959) modeled the flow 
in a laminar viscous trailing vortex for moderate to large Reynolds 
numbers by assuming the axial flow deficit and (unnecessarily) the 
swirl velocity to be small when compared to freestream reference 
conditions. As a consequence of the high Reynolds number and small 
axial perturbation assumptions the governing equations decouple, and 
the analytic solution is the familiar exponential vortex which, by 
virtue of the impulsive flow analogy, is identical to Lamb's (1932) 
for the time development of a viscous vortex. A different class of 
core flow was derived by Hall (1960) as well as Ludwieg (1962) for the 
conical incompressible Euler equations. The chief feature of these 
analytic solutions are radial logarithmic singularities as contrasted 
to the potential (1/r) form. This inviscid solution was extended by 
Stewartson and Hall (1963) to formally account for viscous effects by 
asymptotic means and the resultant solutions exhibit an asymptotic 
viscous subcore. 

Brown (1965) extended Hall's Euler vortex solution to 
compressible conditions. Brown showed that the adiabatic conical 
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governing equations imply the flow to be isentropic and that momentum 
and continuity equations contain the energy integral 

1 q 2 + _JL_E = const. (2.1) 

2 Y-l P 

which, for isentropic flow, will yield the compressible Bernoulli 
equation. The governing equations indicate that the velocity and 
vorticity are locally aligned which is consistant with Crocco's vortex 
theorem for steady inviscid flow: 

Vh Q - TVS = 0 x t. ( 2 * 2 ) 

The chief feature of Brown's solution is that compressibility removes 
the radial singularities at the vortex axis but does so by achieving 
vacuum conditions there. 

Nonaxi symmetric effects for conical flows were studied by Mangier 
and Weber (1966) by modeling a spiral sheet of vorticity in an other- 
wise potential flow. An asymptotic solution was derived, valid for a 
slender core, with the noteworthy result that the lead terms of the 
solution were axi symmetric and identical to Hall's (1960) solution. 
This solution was later extended by Brown and Mangier (1967) to 
account for compressibility yielding a similar result, when compared 
to the axisymmetric solution of Brown (1965). A more general treat- 
ment of the asymptotic structure of spiral vortex sheets including 
nonconi cal effects has been given by Guiraud and Zeytounian (1977). 
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A method valid for both nonconical and conical external flows was 
developed by Hall (1967). The flow is assumed to be steady, axi sym- 
metric, and imcompressible. With the additional assumption of high 
Reynolds number a boundary layer-like approximation valid for a 
slender core is made and the resultant equations, referred to as 
quasicylindrical Navier-Stokes equations, are parabolic in the axial 
direction, i.e., down the core. The equations are solved numerically 
by a finite difference technique and the resultant solutions are shown 
to exhibit many of the general features of vortex core flow for both 
generic wake and leading edge vortex flows. These equations were also 
solved by Bossel (1972) by a method of weighted residuals for wake- 
like initial conditions in a study directed at vortex breakdown 
phenomena. The slenderness assumption was removed by Grabowski and 
Berger (1976) who solved the resultant elliptic system of equations 
with downstream Neumann boundary conditions by an ADI technique for 
initial and edge boundary conditions indicative of wake-like 
vorticies. Their solutions exhibit some features of the vortex 
breakdown flow field which are at least qualitatively similar to 
experiment. 

2.3 Leading Edge Vortex Theory 

In contrast to vortex core formulations, leading edge vortex 
formulations, which are generally directed at wing applications, tend 
to emphasize global flow effects while incorporating very approximate 
representations of the vortex core region. Early work in this area 
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was mostly for conical flows and focused on developing appropriate 

models for the leading edge vortex sheet and vortex core. Contri- 

butors to this research include Legender (1953), Brown and Michael 
(1955), Mangier and Smith (1959), and Smith (1966). Smith's resultant 
model incorporated a segmented free vortex sheet (alligned with the 
local flow and locally force free) terminated with a "cut" or fed 

sheet to a line vortex representati on of the vortex core. By imposing 

that the combination of the fed sheet and vortex core be force free in 
conjunction with no flow boundary conditions for the wing and a 
leading-edge Kutta condition, the problem is solved numerically in an 
iterative fashion. The solutions achieved are qualitatively correct 
but suffer quantitatively due to the assumption of conical flow. 

Several models have been developed for computing the subject flow 
in three dimensions. One approach is to represent the leading edge 
vortex by a system of discrete vortex filaments whose strength and 
location is determined iteratively so that they are force free. 
Contributors to this formulation include Rebach (1973), Mook and 
Maddox (1974), and Kandil, et al . (1974). This formulation has 
generally been coupled with a vortex lattice representation of the 
wing or, in the case of Mehrotra and Lan (1978), a quasi vortex 
lattice. The filament approach generally yields good estimates of 
force-moment properties but published correlations for wing pressure 
distributions are sparce and can be less than satisfactory. An 
alternate approach is to represent the leading edge vortex (as well as 
the wing and near wakes) by higher order panel techniques. The vortex 
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sheets satisfy the usual boundary conditions and the resultant 
formulation takes advantage of the higher order panel methodology to 
provide reasonable estimates of wing pressures as well as resultant 
force/moment properties. This formulation has primarily been 
developed by Brune, et al. (1975) as well as Johnson, et al. (1980) 
and has been applied to a fairly broad class of simple wing shapes. 

The most recent approach to computing leading-edge vortex flows 
is based on solving the three-dimensional Euler equations. The calcu- 
lations are usually founded on a finite volume formulation coupled 
with a fourth-order Runge-Kutta time-stepping scheme as put forth by 
Jameson, et al . (1981). Rotational effects are inherently represented 
in the governing equations, thereby offering the distinct advantage of 
implicitly capturing regions of vortical flow as contrasted to the 
previously discussed methods which (a priori) explicitly model these 
regions. Calculations by Hitzel and Schmidt (1983), as well as by Raj 
(1984), have shown encouraging correlations with experiment; in 
addition, Hoeijmakers and Rizzi (1984) have recently demonstrated an 
encouraging correlation between Euler and "vortex-fitted potential", 
solutions for a 70-degree delta wing at 20 degrees angle of attack. 
However, some controversy regarding the application of these Euler 
methods to the calculation of leading-edge vortex flows has arisen. 

For wings with blunt leading edges, Newsome (1985) recently 
demonstrated for supersonic flow that the Euler solutions are 
extremely grid dependent due to the nature of the artificial damping 
required for numerical stability; this result will hold true for 
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subsonic and transonic applications as well. Of possibly greater 
significance are the occurrence of extensive loses in total conditions 
(entropy production) in regions coincident with the vortical flow, as 
shown by Krause, et al . (1983b) for the sharp leading-edge case. For 
these shockless flows, there is no physically realistic mechanism for 
the entropy production; its source is numerical error. Further 
research into the nature and consequences of these effects will be 
required for a better understanding of leading-edge vortex solutions 
by Euler formulations. 

Additional discussion of computational vortex flow methods may be 
found in the recent and comprehensive review of Hoeijmakers (1983). 

2.4 Vortex Breakdown Theory 

Vortex breakdown analysis can be accomplished either in the 
context of a stability theory or in terms of flow'field analysis. 

Even for the simplified case of an isolated vortex in a tube, vortex 
breakdown (or bursting) has proven to be a very elusive phenomena to 
predict. Although some breakdown theories have enjoyed wider 
acceptance than others, a fully satisfactory and generally applicable 
explanation of this phenomena has yet to be developed. A review of 
many theories directed at this phenomena has recently been performed 
by Leibovich (1983); additional contributions to the problem of vortex 
breakdown can also be found in A6ARD CP-342 (1983). Accurate 
computation of the vortex breakdown flow field has also been hindered, 
chiefly by the turbulent, unsteady nature of this flow as well as by 
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its inherent three dimensionality. Even so, computations based on 
various subclasses of the full governing equations have shown some 
success in capturing various features of this flow field for an 
isolated vortex. 

Ludwieg (1962) derived a stability criterion, often referred to 
as hydrodynamic, by analyzing annular flow with a small disturbance 
theory for spiral perturbations. The stability criterion 

(1 - v* ) (1 - v* 2 ) - (5/3 - v*) w* 2 > 0. (2.3) 

emphasizes the strongly destabilizing effect of radial gradients in 
the axial flow, figure 2.6. Application of this criteria to the 
inviscid, analytic and conical solution of Hall (1960) predicted 
instability when the helix angle tangent, <t>, exceeded 1.16. This 
critical value occurred at the edge of the vortex and was roughly 
confirmed experimentally by Hummel (1965). Ludwieg's criterion 
predicts that the flow will always become unstable to nonaxi symmetric 
modes prior to axi symmetric modes so long as w* f 0. 

The concept of a critical helix angle, <J> C , has been addressed by 

several authors. Bossel (1968) demonstrated, for a simplified vortex 

model with rigid rotation, that small edge disturbances would amplify 

in the core for d> > /2. For a more elaborate model, Bossel (1969) 
c 

found <)> c > 1.92 corresponded to the occurrence of a breakdown bubble 
on the axis of the vortex core. Under rather general conditions. Hall 
(1966) showed, for the inviscid quasicylindrical equations, that 
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gradients in the edge axial flow can be amplified at the vortex core 
axis for a retarding flow. This relationship can be be expressed as 


dw, 

0 

dz~ 


dw‘ 

dF 


- 4 


J 

0 



(2.4) 


For the retarded flow, the integrand is positive. Critical values of 
<t>, therefore, not only promote instability of the edge flow, but also 
greatly promote deceleration of the inner core flow, primarily at the 
axis. A more general treatment of this feature has recently been 
given by Krause (1983) including viscous effects; the amplifying 
feature is retained although modulated by viscosity. 

An alternate interpretation of vortex breakdown has been put 
forth by Benjamin (1962) by drawing an analogy to the hydraulic jump 
from a supercritical state of the flow to a subcritical state. 

Drawing upon variational methods, Benjamin shows that a finite 
transition between the conjugate flows is possible for inviscid 
parallel flow in a cylinder. The supercritical flow corresponds to 
the concentrated vortex and cannot support standing waves whereas the 
opposite is true for the subcritical case. By this theory, vortex 
breakdown corresponds to a transition from the supercritical to the 
subcritical flows. 

Finally, the quasicylindrical Navier-Stokes formulation of Hall 
(1967) can provide an indication of vortex breakdown in a manner 
analogous to the boundary layer equations. As breakdown is 
approached, the flow undergoes rapid axial deceleration and some of 
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the assumptions of the theory are no longer valid. Hall (1967b) draws 
an analogy between this case and the calculation of boundary layer 
flows for conditions resulting in separation. Vortex breakdown, 
therefore, is related to a failure of the quasicylindrical approxi- 
mation and comparisons by Hall (1967b) to the experimental results of 
Kirkpatrick (1965) for a tube vortex show reasonable correlation. The 
predicted breakdown location was several core diameters downstream 
from the measured location consistent with the parabolic nature of the 
theory. Continuation of the calculation beyond the breakdown point is 
impossible by Hall's method. This would require either an inverse 
boundary layer approach or some treatment of elliptic effects such as 
in Grabowski and Berger (1976), or Bossel (1969). 

Additional review material regarding vortex breakdown has been 
given by Hall (1972) and Liebovich (1978). 
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3.0 INITIAL PLANE (CONICAL) CORE FORMULATION 
For the initial plane of the inner vortex core model the outer flow 
is assumed to be representable in some fashion by conical flow 
notions. Matching of this conical solution to the nonconi cal outer flow 
will be addressed in section 6. 

3.1 Intermediate Inviscid Solution 

When the flow is considered to be steady, incompressible, inviscid, 
and axisymmetric the nondimensional governing equations are, without 
approximation 


■kru) + w = 0 

(3.1a) 

r ,r ,z 


2 


uu ,r + wu ,z ' f- - -P,r 

(3.1b) 

u(rv) r + w(rv) , = 0 

9 ' 9 *- 

(3.1c) 

uw, r + ww >z = - P>z 

(3. Id) 


Transforming the independent variables from the cylindrical (r, z) 
system to the spherical-like (t, s) system yields 

-t(w - f) +2|= —"^ (w +'tu) s 

x 1 + t 


(3.2a) 
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-t(w - f)^,t ‘ F - + p »t 


1 + t‘ 


[(w + tu) u >s + tp }S ] 


-t(w - ^-)v t + ^ = — — 
t ’ r 1 + t 


^ = ~ S '~ 2 t(w + tu) v >s ] 


-t(w - -r) w t - tp t = 5 

t 7 »t ,t i + t i 


[(w + tu) w }S + P jS ] 


By further assuming conical flow all s-derivatives are zero and the 
t-derivatives become total derivatives: 

_ t d_ (w u } + 2u = 0 (: 

1 dt ' F t 


t / w ii) du _ v_ = 

_t (w F dt t dt 


* Ux dv (uv) 

-t < w “ F dt t 


, U\ dw _ dp 
< w - t } dt " dt 


Hall (1961) showed that the resultant system of equations may be 
integrated in closed form. The boundary conditions appropriate to the 


in viscid vortex core flow are 


{u = 0; t = 0} 


(v = v , w = w , p - P e ; t = t e } 


(3.4a) 


(3.4b) 
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and Hall's solution, neglecting effects of 0(t 2 ) is 

,, -a, 

2 O 1 t (3.5a) 

W„ 2 

e 

o 9 . 1/2 

S-- [»e - 4 ,n (3 - 6b) 

e e 


H_ - 1 - ^ in <|-) 
e e 


P - P 


w. 


e _ a 2 ln ,t , 12 .2 ,t v 

- = <t> e In (— ) - j a e ln 't~' 


where 




? 1/2 

a g = (1 + 2<t> e ) - 1 


or, equivalently. 


O a p 

^e = a e {1 + Y~ ] 


(3.5c) 


(3.5d) 


(3.6a) 


(3.6b) 


(3.6c) 


The same solution was derived independently by Ludwieg (1962). Hall 
indicates that equations (3.3) can be solved exactly, without the 
slender core (t 2 « 1) assumption, but that the exact solution is 
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algebraically complicated and does not alter the character of the 
approximate solution. The exact solution is presented in Appendix 10. A 
and compared to the slender core solution in figure 3.1. The conse- 
quences of the slenderness assumption are, in general, small; near the 
axis, the nonslender terms are negligible (t/t e << 1) while, near the 
edge, both the slender and nonslender solutions approach the same 
boundary condition values except for the radial flow which has no edge 
boundary condition, but which is small as compared to the other 
velocities. For practical applications, t e < 0.1 and the slender core 
assumption constitutes a very good approximation. 

Several aspects of this solution are noteworthy. Perhaps its most 
prominent feature is the logarithmic singularity in the swirl and axial 
velocities as well as in the pressure. For this distributed vorticity 
solution the singularity is much weaker than for the potential vortex. 
Moreover, the solution never approaches a potential form. 

Further insights are gained by examining the vorticity of this 
solution. From Appendix 10. B, the solution is seen to be highly 
rotational; the axial component of the vorticity is dominant. As 
indicated in section 1.0, the vorticity and velocity vectors for this 
incompressible Euler vortex are locally aligned so that 


0 X t = 0 (3.7) 

A simple calculation, though algebraically cumbersome, can confirm this 
result; vorticity contours are also streamlines. As a consequence, the 
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pressure field of this highly rotational flow is governed by the 
Bernoulli equation and total pressure losses are identically zero. 
These aspects are important to the matching of this flow to the outer 
formulation. 


3.2 Inner Viscous Subcore 

The governing equations for a steady, incompressible, and 
a xi symmetric flow including the effects of viscosity for laminar flow 
a re 



+ w 


,z 



uu 


,r 


+ wu 


,z 


v_ 

r 


+ P 


>r 




zz 



(3.8a) 


(3.8b) 


u(rv) >r + w(rv) >z 

= i*n [ (rv) ,rr ' F< rv >,r + < rv) , Z ;] 
uw r + ww , + p , 


= i— (w + — w + w „) 
Rn v ,rr r ,r ,zz' 


(3.8c) 


(3.8d) 


Evaluating these equations for the slender core conical solution 
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(equations (3.5)) yields 


w ' ^ (H_ + i, 

“e 2 


2z 


zi(Ve!) 

Rn ' _2 ' 


(3.9a) 


,,2 2 
w q 
e e 

2z 


{0} = 


“°2 + ( 1+t2 ) 
2Rn r L W e * 


«e t" + 2 (5- 


(3.9b) 


w p a p 

e e 


{0} = 


i W Q ® Q 

— ( e - ) 

Rn [ z 2 ’ 


(3.9c) 


and continuity is identically satisfied. The nonzero left-hand side of 
equation (3.9a) arises due to the slender core approximation; the 
solution of Appendix A integrates these terms identically. The viscous 
terms for the radial and axial momentum equations arise solely from 
corewise diffusion and will be small for z » Rn" 1 / 2 . The swirl 
momentum equation alone imparts a scale for the viscous subcore which, 
following Stewartson and Hall (1963), is given by 

r = 0 [z(w e azRn)" 1/2 ] (3.10) 

Stewartson and Hall (1963) proceed to construct an inner solution 
for the viscous vortex which will be asymptotic to the outer inviscid 
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solution, equations (3.5). They follow a standard approach from 
singular perturbation theory, the basic steps of which are outlined 
below. By comparing the size of the viscous and inviscid terms, a 
boundary-layer like approximation to equations (3.8) is made by taking 
the limit 

[w e a e zRn •> ® ; (^(WgOgZRn) 1 ^ 2 fixed] (3.11) 

to yield 


•p(ru.) }r + w. >z = 0 (3.12a) 


v 2 



li- p. 

r H i,r 


(3.12b) 

u i ( rv i ) ,r + w i ^ rv i > ,z = 

■R?rE < rv i > , rr 

(3.12c) 

u i w i ,r + w i w i ,z = " p i ,z 

+ E7 ^ w i ,rr + 7 w i ,r^ 

( 3 . 1 2d ) 


where the subscript i denotes an inner solution. The boundary 
conditions appropriate to the inner flow are 

{u.j = v.j = w.^ = 0; r = 0} (3.13a) 

{(v i ,w i ,p i ) T -► (v,w,p) T ; ( 7 ) (w e a e zRn ) 1/2 -► «} (3.13b) 

Equations (3.12) govern the inner viscous flow. In the limit Rn + ® 



24 


equations (3.12) revert to the outer equations (3.1) but with the 
slender core approximation, t^ << 1; the boundary layer like approxima- 
tion (3.11) is equivalent to the slender core approximation in this 
limit. 

Next a new set of independent variables appropriate to the inner 
flow are set down: 


X = In [^(w e a e zRn) 1/2 ] 

n = Y(w e <x e xzRn) 1/2 



x 


i/2 e x 


The outer flow expressed in these inner variables is 


u_ 

w 


" a e r D >.-1/2 
— (VgXzRn) n 


V_ 1/2 [x + ll!UL + I(. in „ . 1 + i-n 

W e A ’ L 2x X C a p 


1/2 


« e x 1/2 [1 + (- In r, + 7 + ^) + •••] 


1 . 1 


(3.14a) 

(3.14b) 

(3.14c) 

(3.15a) 

(3.15b) 

(3.15c) 


w_ 

w 


„ x [1 + 1$JL+ <1- in ,, +i-)] 
e* L Zx X a e 


( 3 . 1 5d ) 
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P - P £ 


w 


_ 2 2 

-i— {1 + + - (-2 In n + 1 + — ) 

2 1 X X a e 


+ lni + Tn_X(. ln , + l + ^) 

4% X e 


+ ^2 [In 2 ti - (1 + —) 1n t)]} 
X e 


(3.15e) 


The expansion for the swirl velocity is valid for 


i!^ui ( . 1ni1 + i +t) l< 1 


(3.16a) 


or 


( l/Z+l/a ) 1/2 (1/2+1/a ) m 

> % ' e x < n < e x e * (3.16b) 


Equations (3.15) are exact. They represent the outer limit of the inner 
flow and suggest the following series expansion for the inner flow: 


p- = a e ( w e a e XZ Rn )" 1/2 Tl A (x»T)) 

e 


(3.17a) 


p- = a e x 1/2 B(x»Ti) 


(3.17b) 
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w i 

= a e xC(x»il) 


(3.17c) 



Y~ agX 2 °(x»Tl) 


with 


(3.17d) 


A(x»t i) = j + A Q (n) + ^ Aj(ti) + + “• Ag (ri) + ... (3.18a) 

B(x» T l) =! B 0 (n) + B^(ti) + '^ B 2^ t1 ^ + *** (3.18b) 


C(x»h) 


In 


2x 


- C. 
X 1 


(n) + 


In 


Co(ti) 


+ ^2 C 3 ( ti) + .. . 
X 


D(x.n) - i + * i DjW * + D 2<”> 

a. A. ^ X ^ 


(3.18c) 


+ hr D^ri) + ••• (3.18d) 

X 

At this point the character of the inner solution and its relationship 
to the outer solution are clear. Stewartson and Hall (1963) proceed to 
cast the boundary conditions (3.13) in terms of the series representa- 
tion (3.18) and to express the governing equations (3.12) for the inner 
flow in terms of the inner variables. The resultant system of equations 
is algebraically complex, but upon equating terms of like order of 
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magnitude, coupled systems of ordinary differential equations may be set 
down and solved in descending order. Initiation of computation is 
expedited by the governing equation for Bq. The solution for Bq 
satifying all boundary conditions can be expressed analytically in terms 
of a confluent hypergeometric function. In addition, Stewartson and 
Hall (1963) cleverly disect the functions of equations (3.18) into 
subsidiary functions which are either universal functions of ti or 
functions of the edge boundary condition term a g . Evaluation of the 
outer limit for the universal q functions is performed numerically and 
the (universal) a e functions are then established once and for all. 

In practice, this allows the inner equations to be uniformly integrated 
in the proper order from the centerline to the edge without iteration. 
Additional details may be found in Appendix C. 

For the present application these equations are integrated with a 
fourth order Runge-Kutta scheme; some typical results are presented in 
figure 3.2. Here the asymptotic nature of the solution is evident as 
well as the consequences of the outer expansion of the swirl velocity, 
equation (3.15c). Unfortunately this approximation renders the 
convergence of the entire formulation to be algebraic as opposed to 
asymptotic, and additional terms in the outer expansion as well as their 
counterparts for the inner expansion would have to be included for 
increased accuracy. 

Because the inner solution is asymptotic the choice of the "edge" 
has no appreciable effect on the resultant solutions so long as 
condition (3.16) is met and the outer and inner solutions agree at the 
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"edge." A comparison of solutions demonstrating the effects of the edge 
location is presented in figure 3.3. For this comparison, the values of 
the t 0 = 0.10 solution at t/t g = 0.05 are imposed as boundary conditions 
for the t e = 0.05 solution and the results are presented as normalized 
by the edge quantities of the t e = 0.10 solution. Small differences 
arise chiefly due to the outer swirl expansion and the definition of 
X. Matching is part of the solution. 

The edge helix angle and the Reynolds number have strong effects on 
the vortex core solution properties. As shown in figure 3.4, an 
increase in edge helix angle, <|> e , greatly increased the magnitudes of 
the velocity components as well as the magnitude of the static 
pressure. An increase in Reynolds number decreases the extent of the 
viscous subcore, much as would be expected (figure 3.5). Reynolds 
number has little effect on the radial flow. However, increases in 
Reynolds number results in increases in the magnitudes of the other 
velocity components as well as the magnitude of the static pressure near 
the axis. Because the radial extent of the inner vortex representation 
(not the viscous subcore) is given by the outer vortex representation 
(section 6), a minimum Reynolds number for which the theory is 
applicable may be derived. This limitation is consistent with the 
assumptions of the outer vortex formulation, namely that the flow is 
inviscid and that free vorticity can be represented by infinitely thin 


free vortex sheets. 
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4.0 N0NC0NICAL INNER FORMULATION 
4.1 Governing Equations 

As in the case of the viscous subcore, the flow is considered to be 
steady, incompressible, axisymmetric, viscous, and laminar, equations 
(3.8). A boundary-layer like stretching is introduced such that 



(4.1a) 

U - >*f u 

(4.1b) 

1 < f « Rn 1/2 

(4.1c) 


and R and U are presumed to be of unit order of magnitude. The 
governing equations (3.8) then become 


h m \R*\z - 0 


fer< uu ,R + wU ,z> - r- - p ,R 


(4.2a) 



,RR 





U(Rv)_ r + v.(Rv)_ z =^[(Rv) >rr 



(4.2b) 


(4.2c) 
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Uw >R + ww >z 




+ 




(4. 2d) 


In the high Reynolds number limit with the constraint of equation (4.1c) 
these equations collapse to the quasicyl indrical formulation of Hall 
(1967): 


ir (RU) ,R + V “ 0 


(4.3a) 

v 2 . „ 

R “ P ,R 


(4.3b) 

U(Rv) #r + w(Rv) >z = (Rv) jRR 

• f (Rv) ,r 

(4.3c) 

Uw,R + ww >z = -P jZ + w >RR + 

1 W ,R 

(4.3d) 


In practice the quasicyl indrical scaling (f = 1) can result in 
values of R and U greater than unit order of magnitude, and the 
scale factor f is useful in analyzing the consequences of these 
effects with equations (4.2). As long as scaling is appropriate 
(equation (4.1c)), it is observed that corewise diffusion effects are 
uniformly negligible among the viscous terms. The remaining viscous 
terms in the swirl and axial momentum equations become less significant 
for f > 1 (much as would be expected); this effect is included in 
equations (4.3c) and (4.3d). The radial momentum equation viscous 
effects are wholly negligible, while the bracketed inertia terms are 
negligible in terms of equations (4.1c). However, these are the 
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"largest" terms of those excluded in the quasi cylindrical approximation, 
and their chief effect is to lessen the magnitude of the radial pressure 
gradient. Nonetheless, the quasicylindrical equations (4.3) are 
reasonable for high Reynolds number flows and are used for this study. 

It should be noted that the quasicylindrical equations which were 
derived with a simple boundary layer stretching notion are identical to 
the initial plane inner equations (3.12) which were derived more 
formally from a boundary layer limit viewpoint. Consequently, the 
equations used to advance the solution in space are consistent with the 
equations used to generate the initial plane solution. 

Following Hall (1967) the equations (4.3) are cast in terms of more 
convenient variables 


K ■ X5v 


. w 
x = w 


to yield 


(m 


.5 


-*,z + S 2 


9 £ 


(4.4a) 

(4.4b) 


(4.4c) 


( 4 . 4d ) 


(4.5a) 


•« xV 


(4.5b) 



I 
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(4.5c) 


ww - 5 - w p + (H - y~ £,w - 

,z x re 


— ) w = — 

2 . w ,e 


~9 ~ P 


3 .2 


x e 


,z 


( 4 . 5d ) 


with 



(4.5e) 


In equation (4.4d) lambda is defined as the ratio of the edge radius at 
any longitudinal station, z m , to the edge radius at the initial 
plane, z^. 

The equations (4.3) or equivalently (4.5) are parabolic in the 
axial direction and therefore can be marched downstream so long as 
appropriate initial and boundary conditions are provided. Although the 
complete flow field in the initial plane is available (section 3) for 
the wing case, only the bounding radius as well as the radial distribu- 
tions of the circulation parameter, K, and the axial flow, w, are 
required to initiate the marching procedure. The boundary conditions 
for the present problem are 

{u = k = w ^ = 0; h; = 0 , z > (4.6a) 

{k = k e , w = w e , p = p e ; l = 5 e , z > z x > (4.6b) 


X(z) given. 


(4.6c) 
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Other edge boundary conditions may be constructed for stream surfaces, 
solid surfaces, or alternate choices of known and unknown edge param- 
eters as discussed by Hall (1967). These alternate conditions will only 
be addressed for appropriate validation cases presented later in this 
section. 

4,2 Numerical Solution Technique 

The solution is advanced in space by replacing the differential 
equations (4.5) with an implicit set of difference equations for each 
computational sub-domain 

f(5.z)|0 < ? < z ra < z < z mtl l 


Second order accurate central differences are used for a grid with 
variable Az m and a simple stretched radial spacing 


A?L 


AC, 


n-1 


(4.7) 


where x is a constant greater than or equal to 1. This permits 
improved radial resolution near the core centerline. Nonlinearities are 
treated with full Newton linearization. The continuity, radial, 
angular, and axial momentum equations are written about points (m + 1/2, 
n - 1/2), (m + 1, n + 1/2), (m + 1/2, n), and (m + 1/2, n), respec- 
tively. Two solution techniques were investigated, the iterative 
decoupled formulation of Hall (1967) and an iterative coupled formula- 
tion based on the Davis modified tridiagonal algorithm as reported by 
Blottner (1979). For the coupled formulation, iteration was required 
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because of the boundary conditions for the continuity and radial 
momentum equations; the continuity, angular momentum, and axial momentum 
equations were coupled while the radial momentum equation was lagged. 

The decoupled formulation was found to converge rather quickly (seven 
iterations, on average) and as such the usual advantages of a coupled 
formulation were offset by the increased level of work per iteration. 

In addition, the run time for a typical core computation is small (* 3 
seconds on CY175) particularly with regard to the outer flow calcula- 
tion. Therefore, the decoupled formulation was used to advance the 
solution. 

Convergence was judged by computing normalized L 2 residuals 
across the entire core for the four computational plane dependent 
variables, H, p, K, and w. The solution at each corewise station 
m was considered converged when all four residuals were of order of 
magnitude 10‘ 4 or less. Additional details may be found in Appendix D. 

4.3 Isolated Core Computations 

Computations were performed for several isolated core cases. The 
purpose of the computations is to validate the present formulation with 
comparisons to published numerical results and exact solutions as well 
as to highlight general features of the computed flow. 

The first case considered was published by Hall (1967) and is 
appropriate to a trailing vortex. The initial conditions at z ^ = 0.25 


are 
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K = 7 (1 - e “ 5 ) 


(4.8a) 


1 -r 

w = 1 - j e 


The boundary conditions at ?, p = 6 are 


(4.8b) 


K " 7 
w = 1 
p = const 
K = 1 


(4.9a) 

(4.9b) 

(4.9c) 

(4.9d) 


The results (figure 4.1) are identical to Hall's (1967) and illustrate a 
diffusion dominated vortex flow with an axial velocity deficit. A 
similar result is shown in figure 4.2 for the same initial and boundary 
conditions but with an initial axial flow excess at z = 0.25: 




(4.10) 


For these solutions, the longitudinal decrease in the axial velocity 
excess (figure 4.2b) exceeded in magnitude the longitudinal increase in 
the axial velocity deficit (figure 4.1b). 

The second case considered, also published by Hall (1967) is 
appropriate to a generic leading edge vortex. Initial conditions are 
derived from the tabulated solutions of Stewartson and Hall (1963) which 
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neglect some of the higher order terms. The edge conditions for 
Zj = 50 are: 


o 

t— 1 

C\J 

11 

0) 

> 

(4.11a) 

w e = 392 

(4.11b) 

£ e = 9.6295 

The boundary conditions are of two types. The first 
a conical "edge" with constant velocity and pressure and, 
increasing circulation: 

(4.11c) 
portion models 
hence. 

50 < z < 100 


v e = 210 

(4.12a) 

• w e = 392 

(4.12b) 

p e = const 

(4.12c) 

= 50 

(4.1 2d) 


The second portion models a stream surface of unspecified shape but with 
a specified pressure distribution, taken in this case to illustrate the 
effects of an adverse pressure gradient. The stream surface will have 
fixed circulation and negligible diffusion across it. These conditions 
necessitate a different treatment at the edge of the core: 
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100 < z < 140 


K e (z) = K fi (100) 

(4.13a) 

p e (z) = prescribed 

(4.13b) 

V , xl 

w Q f, e X 
e 

(4.13c) 

„ dw e x' K e d »e 

w e dz 3 2 dz 

(4.13d) 


The edge parameters w e , H e , and X are unknown, and the solution 
process must therefore be extended to the edge itself. This is 
accomplished with an additional outer iteration procedure where the 
shape parameter X is relaxed and the edge conditions of equations 
(4.13) are satisfied. 

The results (figure 4.3) agree well with those obtained by Hall. 

For the conical zone the solution shows axial acceleration even though 
the edge velocities and pressure are held constant and the core is 
expanding. This can be attributed to the longitudually increasing edge 
circulation boundary condition. The radial extent of the viscous 
subcore does not increase appreciably for this case. For the stream 
tube zone the effects of the adverse pressure gradient are to decelerate 
the flow and expand the stream tube, much as would be expected. 
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The final validation case is for a vortex in a tube with break- 
down. This problem was studied theoretically by Hall (1967b) for the 
experimental results of Kirkpatrick (1964). The test was conducted in 
air with a reference Mach number of approximately 0.08 corresponding to 
a unit Reynolds number of approximately 0.5 million per foot. The 
initial swirl was imparted to the flow by guide vanes and the experimen- 
tal profiles of initial swirl and axial flow can be approximated by 
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The test section of the tube included a small contraction to fix the 
location of breakdown, figure 4.4. 

Stream tube boundary conditions are once again applied for the 
edge, taken as the tube radius. This time, however, the geometry of the 
tube is fixed and the edge quantities H e , w e , p e must be solved for 
to satisfy conditions (4.13c) and (4.13d). Following Hall (1967) these 
edge boundary conditions are incorporated directly into the solution 
process without requiring any additional outer iteration schemes. 

As in the previous cases, the results agree well with those 
published by Hall (1967b), figure 4.5. They demonstrate the type of 
breakdown predicted by the quasicylindrical formulation, a rapid 
deceleration of the axial flow in the vicinity of the vortex axis. As 
the value of the nondimensional axial velocity drops below unity, the 
quasicylindrical assumptions become invalid and the computation will 
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eventually fail. Because the equations are parabolic, the failure tends 
to occur downstream of the experimental location. 

A final correlation is presented for the Burgers vortex, an exact 
solution for the axisymmetric incompressible equations for the following 
assumed form: 

u = u(r) 


(4.15a) 

v = v(r) 


(4.15b) 

w = zw(r) = 2az 


(4.15c) 

The solution is 
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Comparisons between the exact and numerical solutions (figure 4.6) 
show good correlation. This comparison stresses one of the neglected 
terms in the radial momentum equation, uu jP , as well as numerical 


errors. 
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5.0 OUTER FORMULATION - FREE VORTEX SHEET 


5.1 Governing Equations and Boundary Conditions 

The free-vortex-sheet (FVS) theory of Johnson et al . (1980) is a 
fully three-dimensional, inviscid, potential flow formulation which 
models steady subsonic flow about wing or wing-body configurations with 
separation-induced leading- and/or side-edge vortex flow. It is assumed 
that the separation line, a viscous flow phenomenon, is known and that 
boundary layer effects are negligible. Adopting a potential flow 
approach, particularly with regard to the free vortex sheets, is a 
common practice which is supported by experimental studies such as 
Maskell (1964) or Peckham (1961) where the basic shape and strength of 
these vortices were demonstrated to be only weakly dependent on Reynolds 
number (i.e., viscous effects) for slender, sharp-edged wings. The flow 
is governed by the Prandtl Glauert equation 
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(5.1) 


accounting for linear compressibility effects. Here $ is the 
perturbation velocity potential. 

A higher order panel method approach is adopted to solve the 
problem (figure 5.1). The panels incorporate either quadratic doublet 
or linear source distributions. Doublet panels are used to model thin 
wings, wakes, and vortex networks while source panels are used to model 
thickness effects. Thick wings may be modeled with surface source 
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networks and internal mean plane doublet networks. However, most wings 
of interest are sufficiently thin to render the thin wing assumption a 
good approximation which reduces computational expenses. 

The method differs from conventional higher order panel methods by 
virtue of the boundary conditions which must be satisfied to model the 
vortex flow case. The wing must satisfy the usual zero mass flux 
boundary condition locally ("no flow"), but with the Kutta condition 
imposed along network edges with either vortex separation or wake 
flow. The Kutta condition is implemented by matching not only doublet 
strength, but also the component of the vorticity along the appropriate 
edges. This approach, as originally proposed by Rubbert of Boeing, in 
conjunction with the vortex sheet updating procedure, results in all the 
usual flow phenomena at the wing edges associated with the Kutta 
condition, i.e., zero pressure jump, finite flow, smooth off-flow, etc. 

Wakes must satisfy the usual zero pressure jump boundary condition 
locally ("no load"), but with the near wake doing so to second order 
accuracy in perturbation velocity quantities. (The trailing wake is 
formulated in the usual sense and persists far downstream.) This higher 
order formulation is necessary due to the large spanwise velocities 
induced on the wake by the primary vortices, thus rendering the linear- 
ized ACp approximation invalid. A consequence of this higher order 
formulation is that the wake vorticity will skew laterally in the pro- 
cess of satisfying the no-load boundary condition and thereby influence 
the Kutta condition at the wing trailing edge. The near wake accounts 
for the effect which was approximated by the simple transverse wake of 
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Johnson et al. (1976). Inclusion of the near wake has been shown by 
Johnson et al . (1980) to significantly improve correlations between 
theoretical and experimental lift coefficients over a range of delta 
wing aspect ratios. 

The free vortex sheet must have no local flow through it and 
simultaneously support no local pressure jump across it. This duality 
in boundary conditions is balanced by the fact that both the geometry of 
the free sheet as well as its singularity strength are unknown. 

However, by virtue of these boundary conditions the resultant integral 
solution equation is nonlinear. Solution of the subject vortex flow via 
the present approach results, therefore, in a boundary value problem for 
which both the location and strength of a portion of the boundary, the 
vortex sheets, is part of the solution. 

As in conical flow (Smith, 1966), the fed vortex sheet represents 
an approximation to the inner wrappings of the free vortex sheet and 
serves the purpose of terminating free sheet rollup as well as condens- 
ing free sheet vorticity into the vortex core. The fed sheet's strength 
is dictated by the free sheet, while its size and position are deter- 
mined by imposing a boundary condition "consistent with those that would 
be applicable to an infinitely rolled-up vortex sheet as well, namely, 
that the total force normal to the core be zero" (Johnson, 1980) at each 
chord station. Such an approach is analogous to the conical flow formu- 
lation. In practice, free sheet rollup is terminated at approximately 
one half of a revolution, near the crest of the vortex. Sensitivity 
studies for both conical flow and the FVS theory have demonstrated the 
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wing pressures to be reasonably converged at this value of free sheet 
rollup; additional quantities of free sheet rollup do not change the 
solution appreciably. 

Because of the nonlinearity associated with the vortex sheets, an 
iterative solution procedure must be employed. The basic approach is to 
initiate singularity values for a fixed, initial geometry by satisfying 
some but not all of the boundary conditions. The full problem is then 
usually solved by a quasi-Newton scheme with controlled step size. 
Results are considered to be converged when the sum of the squares of 
the residuals is less than 10 . 

5.2 Inviscid Computations 

A correlation with the experimental pressure distributions from 
Luckring, et al . (1982) is reproduced in figure 5.2. Also shown are 
conical vortex flow as well as attached flow estimates. Differences 
between the thoeretical and experimental primary vortex suction peak are 
primarily due to secondary separation effects in the experiment. 

Although the data are for turbulent secondary separation, the secondary 
vortex effects are not small as is evident in the experimental results 
for the outer 20 to 25 percent of the local semispan. Diminishing this 
effect experimentally or adding a representation of the secondary 
separation to the theory would improve the correlation. Even without 
this effect, the correlation is reasonable and substantially better than 
that given by conical vortex flow or attached flow theories. In 
addition, correlation between theoretical and experimental vortex 



44 


trajectories is quite good, as shown in figure 5.3, also from Luckring, 
et al. (1982). 

A correlation with the experimental force/moment properties of 
Hummel, et al. (1972) is presented in figure 5.4 for an aspect ratio one 
delta wing. This sharp-edged wing was beveled asymmetrically with the 
upper surface being flat. The mean plane camber surface effectively has 
a trailing edge flap deflected upward six degrees and the correlation 
between theory and experiment, including this effect, is quite good. 

The free-vortex-sheet theory generally gives good estimates of 
inviscid surface pressure distributions as well as integral properties. 
The formulation has been validated for a fairly broad class of generic 
configurations; these include the geometric effects of leading- and 
trai ling-edge sweep, leading- and trai 1 ing-edge camber, and edge 
cropping, as well as the flow effects of angle of attack, angle of 
sideslip, and freestream Mach number (linear compressibility). Appli- 
cations of this theory have been addressed by Luckring, et al. (1982) as 
well as by Lamar and Campbell (1983), Lamar and Luckring (1979), and, 
most recently, Polhamus (1983). The theory was selected to model the 
outer flow because of the range of documented applications for which 
reasonable pressure estimates have been achieved prior to vortex 


breakdown. 
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6.0 MATCHING - INNER TO OUTER SOLUTION 

The inner flow representation of the vortex core requires "edge" 
values of axial flow, circulation, and pressure as well as the region of 
the edge itself. In the overlap region to be sought between the inner 
and outer theories, both represent inviscid vortices, the former 
modeling axi symmetric, continuously distributed vorticity and the latter 
modeling nonaxi symmetric, concentrated vorticity. Matching requires 
that the two formulations accurately model the flow through this overlap 
zone to a consistent order. For conventional boundary layer flows (two- 
dimensional, axi symmetric) the matching process along with the requisite 
viscous-inviscid interaction has been developed to a fairly mature 
science. However, application of these concepts to the current problem 
is not straightforward; matching must be accomplished in space, away 
from any wing surfaces and the differences between the inner and outer 
formulations must be reconciled. 

6.1 Nonaxi symmetric Effects 

Mangier and Weber (1966) studied the steady incompressible flow 
near the center of a rolled up vortex sheet. They considered a conical 
flow with sheets of concentrated vorticity embedded in an otherwise 
potential flow, which is incompressible. The usual boundary conditions 
are to be satisfied; the sheet must be a stream surface and sustain no 
pressure jump across it. These require the flow to be locally aligned 
with the sheet and, with the additional constraint of constant total 
pressure, that the magnitude of the velocity vector be the same across 
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the sheet. The direction of the velocity vectors across the sheet will 
differ. 

The problem is cast in spherical coordinates by first defining a 
dual stream function to integrate the continuity equation and satisfy 
implicitly the stream surface boundary condition. Next a conical 
potential function is defined and by equating potential and stream 
function representations of the velocity components three coupled first 
order partial differential equations are derived. After transforming 
coordinates and introducing new independent variables an alternate 
system of equations is derived which readily allows the pressure jump 
boundary condition to be implemented exactly; there is no small 
perturbation assumption. The resultant system of equations is solved in 
an asymptotic sense, valid for a slender core. 

The resultant velocity field expressed in cylindrical coordinates 
i s 
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where 9,. is the sheet geometry given by 
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and 
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(6.1e) 


The most noteworthy result of the Mangier and Weber solution is that 
"the leading terms of the velocity components for a potential flow with 
the vorticity concentrated along a vortex sheet are the same as for an 
axisymmetric flow with distributed vorticity" as given by Hall (1961) 
and used herein. In addition, these lead terms of the asymptotic expan- 
sion are continuous across the vortex sheet, and the pressure is 
determined from the Bernoulli equation. To be consistent with the 
axisymmetric vortex core formulation these (higher order) nonaxisym- 
metric effects are neglected in the present formulation and an 
appropriate set of axisymmetric conditions are determined from the outer 
free vortex sheet formulation. The Mangier and Weber solution serves in 
the capacity of an intermediate theory to the matching process of the 
inner and outer formulations. 
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6.2 Matching Criteria 

For conventional free vortex sheet computations the vortex sheet is 
explicitly modeled through approximately n/2 radians of rollup from the 
horizontal, roughly to the crest of the vortex, and terminated with a 
construction commonly referred to as the fed vortex sheet. The fed 
sheet approximates the far field effects of the inner wrappings of the 
free sheet and the amount of rollup explicitly modeled has been shown by 
Johnson et al . (1980) to be sufficient for converging solution 
properties such as wing pressures and vortex geometry much as was the 
case for conical flows (Smith, 1966). A typical result taken from 
Johnson et al . (1980) is reproduced as figure 6.1. It must be 
recognized that the flow properties due to the fed sheet are valid only 
in its far field; flow details in the immediate vicinity of the fed 
sheet may be physically fictitious though mathematically correct. 

The matching between the outer free vortex sheet formulation and 
the intermediate vortex core formulation of Hall is determined numeri- 
cally in an asymptotic spirit; the outer limit of the inner flow and the 
inner limit of the outer flow are studied with the intent of determining 
an overlap zone for which both models describe the flow accurately and 
to a consistent order. Results are presented for the conical initial 
plane formulation at a longitudinal station slightly aft of the apex. 

The conditions chosen for this study are a 70° delta wing at 20° angle 
of attack and a freest ream Mach number of zero. 

The inner vortex properties of the outer free vortex sheet formula- 
tion are determined by systematically increasing the amount of sheet 
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rollup modeled beyond the conventional amount. Solutions were obtained 
in increments of additional rollup of tt/2 radians up to a total rollup 
of 3it radians beyond the datum amount. Typical vortex geometries are 
shown in figure 6.2 for several rollup values. The vortex sheet is 
oblate and its geometry was relatively independent of the amount of 
rollup modeled. The oblateness of the vortex is one manifestation of 
nonaxi symmetric effects. With increased rollup, the radial distance 
from the vortex core to the free sheet decreases and, thus, the 
solutions modeling various extents of modeled free-sheet rollup can be 
used to probe the inner properties of the inviscid vortex at discrete 
intervals. The vortex core properties are determined in a coordinate 
system which locally follows the converged vortex core trajectory as 
given by the free vortex sheet. (See figure 1.1.) The radial extent of 
the vortex is taken as the distance perpendicular to the core trajectory 
from the core to the free sheet terminus; it is measured in the fed 
sheet plane. 

Flow field studies were conducted for the various rollup solutions 
between the free sheet and the core. It was found that, away from the 
immediate vicinity of the free vortex sheet, the circulation was 
constant both in the radial and angular direction and that, as would be 
expected, 

k = ^ (6.2) 

where p c is the doublet strength of the vortex core. The axial flow 
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was found to be constant radially and to vary linearly in the angular 
direction, consistent to the feeding nature of the fed sheet. 

The nonaxi symmetric effects are eliminated from the free vortex 
sheet solutions by taking an angular average at a fixed radial extent 
within the free sheet; e.g.. 



/ wde 

V 2lt 


(6.3) 


as suggested by the form of the nonaxi symmetric terms in the Mangier and 
Weber solution. These properties reflect an inner limit of the outer 
formulation for a given angular extent, and equations (6.2) and (6.3) 
are therefore appropriate for providing boundary condition quantities 
for the vortex core model. These boundary conditions may be interpreted 
as maintaining the same axial mass flow and circulation (as well as 
circulation gradient) for the axisymmetric and nonaxi symmetric flows. 

A comparison of vortex trajectories for the free-vortex-sheet 
solutions with two rollup values and the Mangier and Weber solution 
based on boundary condition quantities extracted from the base line 
free-vortex-sheet calculation (A9=0) is also shown in figure 6.2. The 
correlation is reasonable except for the portion of the sheet given 
by it<A9<2Tt. Radial correlations between the free vortex sheet and Hall 
inviscid solutions are presented in figure 6.3 for the axial and swirl 
velocities as well as the longitudinal circulation gradient. The base 
or "edge" solution is taken as the free vortex sheet solution with 
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either the standard amount of rollup (solid line) or an additional 
amount of rollup of %/Z (dashed line). The correlation is quite good 
among the points in the fed sheet plane, but not as good among the 
points it/ 2 out of phase. The difference between these two sets of 
results is viewed in part as a remaining nonaxi symmetric effect most 
probably due to the oblate nature of the free vortex sheet. Additional 
studies using the solutions with % and 3it/2 of additional rollup as the 
"edge" showed the correlation in the base fed sheet plane to be equally 
good while the correlation in the plane % out of phase was still not as 
good. The character of these solutions together with the results of 
Mangier and Weber suggest that an asymptotic approach could account for 
these nonaxi symmetric effects. Because the calculations in the plane of 
the fed sheet with no additional rollup correlated best with the inner 
formulation, they were chosen to provide the appropriate boundary 
conditions. It is to be emphasized that within this plane the overlap 
is independent of the amount of rollup modeled. 

Consistent to the theoretical studies of Hall and Mangier and' 
Weber, the free vortex sheet boundary conditions result in a pressure 
field which is continuous across the vortex sheet and given by the 
Bernoulli equation; there are no total losses. Specification of the 
velocity field therefore implies the pressure field and the requisite 
boundary conditions for the vortex core from the free vortex sheet are 
complete. 

It must be emphasized that the particular boundary conditions 
chosen for this initial formulation are not the only ones available. 
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Additional studies should yeild improved correlations as implied by the 
excellent correlation achieved by Kuchemann and Weber (1965) by fitting 
the Mangier and Weber solution to an experimentally determined 
trajectory. The correlation was not limited to the inner wrappings but 
extended well into the outer portion of the vortex. A more rigorous 
treatment of the subject boundary conditions could also entail quantifi- 
cation of nonaxi symmetric effects, a candidate topic for future studies. 

With the exception of radial flow, the solutions are fully 
matched. As in boundary layer analysis, the edge value of radial flow 
may not be specified but must be accommodated through an appropriate 
viscous/inviscid interaction scheme. Various schemes have been studied 
to adjust the overall flow in such a fashion as to fully match the inner 
and outer representations. Lighthill (1958) studied several approaches 
to account for the viscous zone effects on the outer flow and a fairly 
recent review of viscous-inviscid interaction may be found in AGARD CP- 
291 (1981), in particular, in the papers by Melnik (1981) and Lock 
(1981). For the present formulation these effects are not addressed. 
Consistent to the assumed vortex core model, the radial flow is small 
and, for most cases, will have a correspondingly small effect on the 
overall solution. Calculations by Hoeijmakers et al . (1979) showed that 
the additional entrainment increments as modeled by free line 
source/sinks along the vortex core have a small effect on the free 
vortex sheet geometry while affecting the vortex strength to a larger 
extent. Radial flow effects will become increasingly important as the 
axial flow gradients become large, such as in the vicinity of a burst as 
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modeled by the Hall theory. Here the radial flow effects, properly 
accounted for, would cause the vortex to flair and the longitudinal 
pressure gradients to increase as suggested by Bossel (1968) thus 
promoting the upstream propagation of the burst point. Inclusion of 
radial flow matching by an appropriate viscous-inviscid interaction 
scheme is left as a candidate topic for future studies. 

6.3 Nonconical Effects 

Matching the initial -plane formulation to the inner solution 
requires detailed consideration of the nonconical aspects of three- 
dimensional subsonic flows near the apex of a slender wing. Conical 
flow theories such as those due to Jones (1946) for attached flow or 
Smith (1966) for leading-edge vortex flow are often used subsonically as 
a qualitative datum to judge three-dimensional effects, most commonly 
the trai 1 i ng-edge flow effects affiliated with the Kutta condition. 
However, the subsonic flow near the apex of a delta wing is also 
fundamentally nonconical. Sketches highlighting this nonconical nature 
of the flow are presented in figure 6.4 for attached flow on a delta 
wing. The sketches illustrate trends in pressure coefficinet and are 
indicative of other flow properties as well. From symmetry considera- 
tions the isobars down the wing center line are normal to it and, as 
indicated by Thwaites (1960) among others, this condition results in 
large isobar curvature near the apex region. As a consequence, wing 
properties take on a singular form along a ray in this region. 
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The nature of this singularity has been studied for the semi- 
infinite delta wing by Germain (1955) who showed that the linearized 
apex flow is dominated by eigensolutions of the form 

$ = s 1+e <t>(<M) ; -l<e<0 (6.4) 

where <|> is a perturbation velocity potential, s is the distance from 
the apex, and e is a measure of the singularity strength. Values of 
s have been determined numerically by Rossiter (1970) who also reports 
several series expansion solutions valid for leading-edge sweep angles 
near zero or 90 degrees; e is generally small for conventional delta 
wing sweep values. The apex singularity is weak and similar in form to 
two-dimensional wedge flow (see Van Dyke (1964)). 

The nonconical nature of this solution is manifested as a uniformly 
scaled flow. For example in sphericl coordinates the velocity 
components 

(6.5a) 
(6.5b) 
(6.5c) 

take on the form 



v s = (1 + e) s c *(<M) 


(6.6a) 
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v = s'* . 
<P ><P 


(6.6b) 



1 
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(6.6c) 


For e = 0 the conical form is recovered. However, even for e + 0 
ratios of velocities have no s dependence. 

Flows in the general class of equation (6.4) are weakly singular 
along a ray and take on the form 

2 p 

s e = 1 + £ In s + y ~ In* 1 s + ... (6.7) 

This term contributes the major nonconical effect. Consider, for 
example the flow with continuously distributed vorticity governed by 
equations (3.2). If the velocities take on the form 

$(s,t) = s e $(t) (6.8) 

the subject governing equations (3.2) take on the form 
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These equations exhibit no s dependence but differ from the conical 
equations (3.3) by virtue of the right-hand side being of 0(e) as 
opposed to being equal to zero. Equation (6.9) along with equa- 
tion (6.7) suggest an asymptotic solution with the lead solution given 
by Hall, equations (3.5). A similar result can be demonstrated for the 
conical solution of Mangier and Weber. Closed form solutions to this 
class of flow and not pursued in this report. 

For the present problem, the free-vortex-sheet theory removes the 
strong leading-edge singularity with the leading-edge Kutta condition. 
However, the weak apex singularity remains and its form, as determined 
from numerical experimentation, is similar to the attached flow case. 
Representative velocity boundary condition values are presented in 
figure 6.5. Because of the particular form of the apex singularity, 
matching the three-dimensional outer flow to the conical flow initial 
plane solution can be accomplished in terms of locally scaled vari- 
ables. An example is presented in figure 6.6. Shown on the left 
portion of this figure are radial distributions of inviscid axial flow 
in the vortex core at three-chordwise stations near the apex of a 70° 
delta wing at 20 degrees angle of attack. These were obtained by 
matching the intermediate solution to the free vortex sheet boundary 
conditions as previously described. Shown on the right portion of the 
figure are the same results, each scaled by the value of the edge axial 
flow for the three longitudinal stations. In this locally scaled plane 
there is virtually no dependence on longitudinal position and the 
conical initial plane solution can be implemented without constraint. 
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Numerical studies have verified the solutions to be virtually 
independent of the initial plane location. 

6.4 Downstream Continuation of Boundary Conditions 

As described in section 5, the near wake is modeled 2 1/2 to 3 
vortex diameters at the trailing edge downstream. By this point the 
vortex is essentially fully fed and the corewise gradients of 
circulation and axial flow boundary conditions are nearly zero: 

K e,s “ w e,s “ 0 (6,10) 

Continuation of the free vortex sheet modeling much beyond this distance 
can result in convergence difficulties because the strength of the sheet 
is becoming very small and (unmodeled) trailing edge wake rollup effects 
are becoming more significant. However, as described by Hall (1967b) 
the bursting phenomena represented by this vortex core model can occur 
several vortex diameters downstream of the experimental location due to 
the parabolic nature cf the governing equations; in an analogous fashion 
simple boundary layer methods can separate several boundary layer thick- 
nesses downstream of the correct location. Consequently, it is 
desirable to continue the vortex core calculation downstream of the 
modeled near wake extent. Free vortex sheet quantities in this zone are 
inadequate for these purposes because the modeling is that of potential 
flow trailing wakes and valid only in a far field sense. 

Boundary condition quantities for this zone were obtained from the 
free vortex sheet by means of analytic continuation. By virtue of the 
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diminished free sheet strength in the aft portion of the near wake, the 
radial gradients of the circulation and axial flow near the edge of the 
vortex core are also virtually zero. In this portion of the vortex core 
the flow has become that of a potential vortex, and the radial extent or 
"edge" of the core no longer affects the calculation. The boundary 
condition quantities K e , w e , and R 0 are extended downstream by first 
establishing their values and nominal corewise gradients several sta- 
tions upstream from the trailing wake to minimize any near field errors 
affiliated with the trailing wake itself. Using the vortex core diam- 
eter at this location, the quantities are then extended several 
diameters downstream by maintaining continuity of corewise derivatives 
and stipulating that these derivatives go to zero. From this location 
the calculation are continued downstream by maintaining the zero 
gradient condition. The edge pressure always follows from the Bernoulli 
equation. Numerical studies demonstrated that the extent over which the 
analytic continuation was employed had no appreciable effect on the 
computed core properties. 
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7.0 RESULTS AND DISCUSSION 

Calculations for the composite formulation have been performed for 
delta wings ranging in leading-edge sweep angle from 55 degrees to 85 
degrees over an angle-of-attack range from 5 degrees to 50 degrees. 
Select arrow and diamond wings have also been analyzed. All composite 
calculations have been performed for incompressible conditions at 
Reynolds numbers on the order of one million based on the wing root 
chord. Some results highlighting the character of the composite 
solutions are presented below. 

7.1 General Flow Features 

The core flow properties of a 75-degree delta wing at 15 degrees 
angle of attack are presented in figure 7.1 as a function of radial 
distance from the core axis normalized by the local edge radius for 
various longitudinal stations. The radial velocity (figure 7.1a) 
distributions remain essentially linear and decrease in magnitude down 
the core. In figure 7.1b, the circulation parameter radially increases 
at all stations throughout the core. Near the axis, the quadratic 
region, indicative of solid-body rotation, is very small. At the edge, 
the inner solution matches to the outer solution without necessarily 
achieving zero slope radially; a potential -like vortex flow begins to be 
evidenced towards the edge of the most downstream station shown. As 
part of the solution, this inviscid flow transitions from the highly 
rotational leading-edge vortex structure to the irrotational wake-like 
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vortex structure; this transition conies about as the feeding of the 
vortex core from the free vortex sheet goes to zero. 

These effects are also manifested in the angular velocity (figure 
7.1c) as well as the axial velocity (figure 7 . Id ) . Near the edge of the 
most downstream station shown, the angular velocity begins to evidence a 
potential (1/r) form while the axial flow is becoming constant. 

Although the edge flow is decelerating, the angular and axial velocity 
components near the axis show a longitudinal acceleration over the 
forward portion of the wing followed by a deceleration over the aft 
section of the wing as well as the wake. In addition, the normalized 
radial extent of the viscous subcore decreases downstream (Fig. 7.1c). 

Consistent with the velocity profiles, the static pressure 
coefficient (figure 7.1e) shows an adverse longitudinal pressure 
gradient at a fixed fraction of core radius away from the axis. Both 
proverse and adverse longitudinal pressure gradients are evidenced in 
the vicinity of the axis. The total pressure coefficient distributions 
(figure 7. If) emphasize the confinement of total losses to the viscous 
subcore by this theoretical method. It should be noted, though, that 
the magnitude of the total losses for this flow are less than the total 
losses associated with simple boundary-layer flow because of the flow 
acceleration which occurs in the viscous region of the vortex core. The 
erroneous overshoot in the initial plane of the total pressure 
coefficient is seen to quickly damp out. 

These same results are presented in figure 7.2 for an unsealed 
abscissa. The positive radial flow, figure 7.1a, occurs at the 
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trai ling-edge station as a consequence of comparatively large 
longitudinal retardation of the axial flow in this region. For this 
flow, the stream tubes will exhibit a bulge in the vicinity of the 
trailing edge. In figure 7.2b the circulation distributions are seen to 
roughly overlay one another; the circulation for an inner portion of the 
core is chiefly established by conditions upstream in these solutions. 
Related to this effect, the radial extent of the viscous subcore (figure 
7.2c), though increasing, changes very little for the range of computa- 
tions shown. At the downstream stations, the angular and axial velocity 
components, as well as static and total pressure coefficients, show 
little dependence on axial distance for the outer portion of the core. 
However, the entire core flow over the forward three-quarters of the 
wing is highly three dimensional as is the inner portion of the core 
flow over the entire computed range. 

Longitudinal variations for select flow parameters are presented in 
figure 7.3. The nonconical and, in general, three-dimensional nature of 
the solution is evident. Even though the edge boundary condition 'for 
the axial flow is continually decelerating (figure 7.3b), the centerline 
value accelerates over the forward portion of the wing. This is due to 
the vigorous increases in circulation which the vortex experiences 
longitudinally in this region; the circulation and axial flow are 
closely coupled by the governing equations, chiefly through the 
pressure. Accordingly, the static pressure coefficient exhibits a 
proverse pressure gradient over the forward portion of the wing even 
though the edge condition is uniformly adverse (figure 7.3c). Losses in 
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total pressure accummulate down the length of the wing along the 
centerline and appear to be approaching an asymptote by the last station 
shown (figure 7.3d). 

The effects of angle of attack and leading edge sweep angle are 
summarized in figure 7.4 and 7.5, respectively, at a longitudinal 
distance down the core of one half root chord, z = 0.5. With increasing 
angle of attack the vortex is strengthened and, correspondingly, 
velocities increase through the core while pressures become more 
negative. Angle of attack has a strong effect on the magnitude of the 
velocities, particularly near the vortex core axis. For the largest 
angle of attack shown, this flow is already approaching a condition 
where compressibility effects are no longer negligible and for angles of 
attack approaching those conducive to near field vortex breakdown, 
compressibility effects will become increasingly significant. An 
increase in angle of attack is also seen to increase the radial extent 
of the viscous subcore over which total pressure losses are manifested, 
(figure 7.4d). 

An increase in leading edge sweep angle has the opposite effect of 
an increase in angle of attack; the vortex is weakened and, correspond- 
ingly, velocities decrease through the core while pressures become more 
positive due to this effect. The solutions of figure 7.5 show essenti- 
ally the opposite trends of the solutions of figure 7.4. 



63 


7.2 Experimental Correlation 

Calculations for the 75° delta wing at 15° incidence are compared 
to the experimental results of Earnshaw (1962) in figure 7.6. The test 
was conducted with an aspect-ratio-one delta wing (76° leading-edge 
sweep) at 14.9 degrees incidence with a freestream reference Mach number 
of approximately 0.09 and a Reynolds number of approximately 3 million. 
The measured velocity and pressure profiles through the vortex core were 
obtained with a five tube Conrad probe aligned to the local flow and 
were taken along mutually perpendicular traverses in the wind axes at 
70-percent root chord. These data were also used by Hall (1960) to 
correlate an earlier vortex core theory. This predecessor theory was 
conical, and the intermediate formulation was empirically fit to the 
experimental results for the boundary conditions at the edge of the 
vortex. In addition, the viscous subcore was treated in a simpler 
fashion than the Stewartson and Hall theory. 

Comparisons between the theoretical and experimental velocity 
profiles show reasonable correlations for the outer region of the vortex 
core and for the radial extent of the viscous subcore. The major . 
discrepancy of this correlation is the centerline axial flow, and both 
theory and experiment are probable contributors to this case. Although 
the probe was very small compared to the local wing semispan, its 
diameter was still appreciable compared to the scale of the viscous 
subcore. Apart from perturbing the flow itself, gradients across the 
probe head will also affect the measurements. In addition, the 
relatively small probe size was achieved in part by use of a rather 



64 


large model set at high incidence relative to the test section, thereby 
creating possibly large wall interference effects. Theoretically, the 
major factors affecting the subcore flow are probably the incompressible 
and laminar flow assumptions. Because these flows have a local maximum 
in velocity at or near the vortex axis, they can be locally compressible 
at incompressible reference conditions. The current solution is 
approaching this condition, and further angle-of-attack increases will 
aggravate the incompressible assumption. The compressible flow would 
undergo less acceleration than the incompressible counterpart, an effect 
which has been studied by Brown (1965) for inviscid conical flows. With 
regard to the laminar flow assumption, the implementation of a simple 
eddy viscosity model or, for that matter, performing the computation at 
a lower Reynolds number, will result in reduced centerline axial flow 
due to the correspondingly increased radial diffusion effects. Hall 
(1960) put forth this argument as the chief reason for the similar lack 
of correlation between his conical theory and this data set. 

The correlation between theoretical and experimental pressure 
coefficients (figure 7.6b) is consistent with the velocity corre- 
lations. The centerline static pressure coefficient is more negative 
than the experimental value chiefly because of the increased axial 
flow. However, static pressures of this magnitude at the core axis are 
not unusual for vortex core flow. At a higher angle of attack Earnshaw 
(1962) recorded C p values of approximately -24; for a 65° swept wing 
at 15° incidence Lambourne and Bryer (1960) recorded C p values in the 
vicinity of -13. The theoretical total pressure losses are confined to 
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the viscous subcore, whereas, experimentally, they are spread over the 
majority of the region shown. Further understanding of the differences 
between theoretical and experimental results will require correlation 
studies for additional configurations and flow conditions. 

7.3 Vortex Breakdown Analysis 

As stated in the background and review section, the variety of 
theories directed at vortex stability and, in particular, vortex 
breakdown, have yet to yield satisfactory correlations with experiment, 
even for the isolated vortex case. For the present analysis, two 
criteria of differing types are selected to evaluate the subject flow, 
the bounday layer analogy of Hall (1967b) and the stability theory of 
Ludwieg (1962). In addition, an empirical analysis of the solutions is 
performed in terms of the maximal helix angle tangent, <t>max» a 
parameter to which breakdown has been shown theoretically to be 
sensitive. 

A common attribute of vortex breakdown is the rapid axial 
deceleration of the core centerline flow, which can lead to an axial 
stagnation point. Under these conditions, some of the assumptions 
leading to the quasicyl indrical equations are no longer valid. As 
pointed out by Hall (1967b), this is in close analogy to the calculation 
of boundary-layer flows for conditions resulting in separation and, 
therefore, vortex breakdown could be related to a failure of the quasi- 
cyl indrical approximation. Hall proceeded to quantify this observation 
and included a correlation to the experimental results of Kirkpatrick 
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(1965). The calculations were performed at l/50th the experimental 
Reynolds number as a means to simulate the turbulent experiment in the 
laminar theory by an eddy viscosity concept. Breakdown was predicted 
several core diameters downstream of the experimental location 
consistent with the parabolic nature of the formulation. This 
particular criterion emphasizes flow properties in the viscous subcore 
of the vortex. 

For the present delta wing applications, no breakdown of this type 
was observed. Although the axial core flow is decelerating, its 
magnitude at breakdown conditions, as determined by Uentz and Kohlman 
(1968), is still sufficient to prevent the quasi cylindrical equations 
from failing. The computations were performed without Hall's eddy 
viscosity scaling because the status of the core flow (laminar or 
turbulent) is uncertain and because compressibility will also affect 
this flow by reducing the axial magnitudes. 

In contrast to the boundary-layer analogy of Hall (1967b), 

Ludwieg's hydrodynamic stability criterion is derived from flow 
stability considerations, is inviscid, and tends to be manifested near 
the edge of vortices. An application of this criterion is shown in 
figure 7.7 for a 70° delta wing at 30° angle of attack. Experimentally, 
Wentz and Kohlman (1968) determined <xb - 29° for this wing. In this 
figure, contours of the stability relationship in the vortex core are 
presented; the roughly diagonal edge where the contours terminate 
corresponds to the edge of the inner computational space as given by the 
free-vortex-sheet theory. In the shaded region, the flow is unstable by 
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Ludwieg's criterion. The zone differs in character from the conical 
case in that it does not correspond to maximum values of helix angle. 

It was found for other wings and other angles of attack, but, in 
general, it did not move so far upstream as the trailing edge. Two 
aspects of the present formulation which would inhibit the forward 
propagation of the unstable zone are the parabolic assumption and the 
lack of viscous-inviscid interaction effects. 

The final criterion chosen for analysis of vortex breakdown was the 
occurrence of a critical helix angle in conjunction with a longitudi- 
nally retarding flow, or equivalently, a longitudinal adverse pressure 
gradient. As put forth by Hall (1972), these two effects, along with a 
divergence of stream tubes in the vortex core, commonly occur in 
conjunction with vortex breakdown. The critical helix angle is also a 
consequence of simplified stability analyses, such as those due to 
Ludwieg (1962) or Bossel (1968), as well as more advanced flow field 
analyses given by Bossel (1969, 1972). 

Contours of <{> in the vortex core are presented in figure 7.8 for 
a 70° delta wing at an angle of attack below, essentially at, and above 
at), respectively, as determined experimentally by Wentz and Kohlman 
(1968). Several aspects of these solutions are noteworthy. For a given 
angle of attack, there are two regions of maximum <j>. Both are at the 
edge of the vortex, one near the apex and the other in the vicinity of 
the trailing edge. The structure of the vortex in terms of this 
parameter changes from the wing to the wake; over the wing, the maximum 
4 occurs at the vortex edge, whereas, in the wake, this maximum occurs 
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well within the vortex. Finally, the maximum 4> is seen to exceed 
unity as ab is surpassed, figure 7.8b. 

In terms of quantifying the condition of incipient vortex breakdown 
at the wing trailing edge, the longitudinal regions of maximum $ are 
considered. Longitudinal pressure gradients for this wing are 
qualitatively similar to those shown for the 75° delta wing (figure 
7.3). Near the apex, the flow may achieve critical <)>, but the core 
axis is in a favorable pressure gradient. Near the trailing edge, 
however, the local maxima of <j> occur in conjunction with an adverse 
pressure gradient down the core axis. The combination of these effects 
along with the low total energy of the core down the centerline could be 
conducive to vortex breakdown. 

The various solutions were analyzed for these conditions at the 
trailing edge. The flows were in general decelerating (adverse pressure 
gradient) and exhibited maximum $ at the edge of the vortex; these 
values are presented in figure 7.9 for several of the delta wings 
analyzed. For these wings, the experimentally determined value of ab 
from Wentz and Kohlman (1968) is also indicated. The experimental 
condition of vortex breakdown at the trailing edge of the delta wings 
roughly correlates with a constant theoretical value of He, max* With 
the present formulation, this value is slightly greater than one; for 
reference purposes, the conical value of 1.16 is also shown in figure 

7.9. 

This correlation is shown in figure 7.10 for a more familiar 
parameter space; the theoretical results of this figure are based on 
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<(> c = 1. In addition to the reasonable correlation with delta wing 
leading-edge sweep effects (which are generally strong), the correlation 
with trailing-edge sweep effects (which are generally weak) is also 
shown to be reasonable for 70° arrow and diamond wings. The significant 
aspect of this correlation is that the single concept of a critical 
helix angle in conjunction with the adverse pressure gradient provides 
reasonable estimates of both the strong leading edge and the weak 
trailing-edge effects and does so with one value of <j> c . The particular 
value of <}> c is less significant. It should also be noted that <t> c 
occurs in a region of the vortex core model where three-dimensional 
effects are strong and, in addition, where the chief theoretical concern 
is the matching of inner and outer representations. 

In light of these results, the lack of a physically realistic burst 
in the computation can be attributed to neglected effects, many of which 
affect the subcore flow to a much greater extent than they do the rest 
of the core flow. Virtually all of these, when considered individually, 
are conducive to the occurrence of a physically realistic burst. 
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8.0 RECOMMENDATIONS FOR FUTURE STUDIES 

Several effects currently not included in the leading-edge vortex-, 
core formulation could readily be implemented. The method could be 
extended to account for compressibility in a straightforward fashion. 
This would account for the local effects near the vortex axis, reducing 
the magnitude of the velocities there. Implementing a simple eddy- 
viscosity model to approximate turbulence effects would also result in 
decreased velocity magnitudes in this region. Of more importance, 
though, would be the implementation of a viscous/inviscid interaction 
scheme to fully match the inner and outer flow models. This will be 
most important under conditions conducive to vortex breakdown where 
deceleration of the axial flow in the core must be accompanied by a 
radial expansion of the vortex. This expansion would tend to promote 
the upstream propagation of a burst. 

Additional matching studies should yield improved correlations of 
the vortex sheet trajectories between the inner and outer representa- 
tions of the vortex. In addition, nonaxi symmetric effects may be ‘ 
addressable for the parabolic formulation. 

Application, of the present approach to computing vortex flows is 
not strictly limited to the methods currently in use. For example, 
elliptic effects could be addressed by embedding a form of the Grabowski 
and Berger (1976) model of the vortex core into the free vortex sheet. 

In addition, other vortex flow problems such as forebody vortex effects 
of importance to missile aerodynamics could be addressed in a similar 
fashion with models appropriate to that flow. 
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9.0 CONCLUSIONS 

The flow in the core of a three-dimensional separation-induced, 
leading-edge vortex can be calculated by appropriately matching inner 
and outer representations of the vortex. This procedure is not strictly 
limited to the theories implemented in the present formulation; it is 
applicable for other representations of the inner and outer flows as 
wel 1 . 

The calculated results of the present formulation exhibit many of 
the general features of the subject flow. Some prominent aspects of 
these solutions can be summarized as follows: 

1. At practical Reynolds numbers, the majority of the core flow is 
inviscid, rotational, and weakly singular in the radial direction. 

2. Near the vortex axis, the flow exhibits velocity magnitudes 
several times the freestream reference value along with total pressure 
losses; these losses arise because of modeled viscous effects. 

3. The solutions are highly three-dimensional and begin to 
evidence a transition in the near wake from the leading-edge class of 
vortex flow to a more wake-like vortex flow structure. 

Although the experimental comparison was limited in scope, the 
correlation with measured vortex core properties is encouraging. 

However, extensions to the model may be required for more accurate 
estimates of the viscous subcore flow, particularly under conditions 
conducive to vortex breakdown. 

The experimental condition of vortex breakdown at the trailing edge 
of delta, arrow, and diamond wings was found to correlate closely with 
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the theoretical criterion of a critical helix angle in conjunction with 
an adverse pressure gradient. The calculations also exhibited a zone of 
unstable flow by Ludwieg's criterion, but did not exhibit failure of the 
quasicyl i ndrical approximations. The lack of an axial stagnation point 
or other realistic aspects of the vortex breakdown flow field is related 
to approximations chiefly affecting the viscous subcore. Systematic 
extension of the present formulation should provide additional insights 
to the vortex breakdown phenomena for three-dimensional flows. 
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10.0 APPENDICES 


10. A Exact Solution for Hall 's Vortex 


The governing equations (3.3) are appropriate for steady, inviscid, 
incompressible, and conical flow. The solution procedure of Hall (1961) 
may be used with more complicated integrals to yield the following exact 
solution subject to boundary conditions (3.4): 
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and, equivalently. 
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The relationship of this exact solution to Hall's slender core 
solution (t 2 « 1) is apparent with the following expansions: 
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and the slender core solution is recovered. 
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10. B Vortical Aspects of Hall's Vortex 


For axi symmetric flows 
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Near the axis (t « 1) the vorticity has the form 
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10. C Stewartson and Hall Inner Solution 

Transforming equations (3.12) to the inner dependent and 
independent variables results in the following coupled system of 
ordinary differential equations written in descending order: 


0 ( 1 ) 


nAg + 2A 0 = -nC^ 


n* - " 2 r 2 

D 1 - — B o 


B 0 + *2 " + ri’ B 0 ' *2 B 0 " 0 

V) 


q + (1 » + c i - r 


tlDJ 


0 (15-i) 

X 


r)A^ + 2A ^ = — riC^ 


°2 ■ T B 0 B 1 


B i' + ( ? 11 + B i ' h B i = r ” B o 

n 


c 2 + n + C 2 = J ^2 “ T 


( C . 1 a ) 
(C.lb) 


(C.lc) 


(C.ld) 


(C.2a) 


(C.2b) 


(C.2c) 


(C.2d ) 



78 


nAg + 2Ag = \ riC^ - tiCj 
D 3 ■ T B 0 B 2 

B 2 + ( 7 " + n* B 2 ' T" B 2 = ^ " C 1 B (> 

' 11 

“ 1 A o B o tiA o b o 

C 3 + 4 ^ + ^ C 3 = T ^3 “ T D 1 " k 


(C.3a) 


(C.3b) 


(C.3c) 




tlCJ(C 1 


+ V 


( C . 3d ) 


Here the prime denotes a total derivative with repect to r\. The 
boundary conditions are 
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The treatment of edge boundary conditions is further simplified by 
introducing subsidiary functions: 
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10. D Numerical Aspects of Inner Core Solution 

Equations (4.5) are evaluated in difference form about points 
(m + 1/2, n - 1/2), (m + 1, n + 1/2), (m + 1/2, n), and (m + 1/2, n) 
respectively. The difference approximations are formally second order 
accurate: 
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The superscript ( ) denotes a value from the previous iterate. For the 
centerline the coefficients of equations (D.13) take the form 


AU = 


-8 


A Wl + 


= - AU 


AD 


AZ 


m 


AL = 0 


ws 3 ,111 \ P m+l,n “ p m,n 

AR ’ m; - AU( Vn+l ■ Vn> 


AZ 


(D.23) 

(D.24) 

(D.25) 

(D .26) 
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The difference equations are valid for the decoupled iterative 
solution procedure of Hall (1967). With an estimate of the axial flow 
at station (m + 1), the continuity equation is integrated for the radial 
flow parameter H. Next the swirl momentum equations are solved 
simultaneously for the circluation parameter K with the Thomas 
algorithm. With this result the radial momentum equation is integrated 
for the pressure. Finally the axial momentum equations are solved 
simultaneously for a new estimate of the axial flow. Convergence is 
judged as described in the text. 
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Figure 1.1- Theoretical flow model. Rn » 1. 
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b) Near vortex core 


Figure 1.2- Frequently used parameters. 



Figure 2.1- Vortex flow structure and lift effect - slender wings. 
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Figure 2.2- Near field vortex breakdown. 
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Figure 2.3- Incompressible vortex flow regimes for sharp-edged 
delta wings. 
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Figure 2.4- Vortex breakdown and lift effects 
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Figure 2.5- Vortex breakdown moment effects. a = 65° 
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Figure 2.6- Hydrodynamic stability criterion - Ludwieg. 
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a) radial velocity 

Figure 3.1- Comparison of slender and exact inviscid solutions. 

Rn = 10 6 , <j> e = 1. 
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b) angular velocity 
Figure 3.1- Continued 
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c) axial velocity 
Figure 3.1- Continued 
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d) static pressure 
Figure 3.1- Concluded. 
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a) radial velocity 

Figure 3.2- Inner and outer initial plane solutions. 
Rn = 10 6 , <(> e = 1. 
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b) angular velocity 
Figure 3.2- Continued 
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c) axial velocity 
Figure 3.2- Continued 
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d) static pressure 
Figure 3.2- Concluded. 
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b) angular velocity 
Figure 3.3- Continued 
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c) axial velocity 
Figure 3.3- Continued 
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d) static pressure 
Figure 3.3- Concluded. 
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a) radial velocity 


Figure 3.4- Edge helix angle effect. Rn 
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b) angular velocity 
Figure 3.4- Continued 
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c) axial velocity 
Figure 3.4- Continued 
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d) static pressure 
Figure 3.4- Concluded. 
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c) axial velocity 
Figure 4.6- Continued 


140 



d) static pressure 
Figure 4.6- Continued 
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e) Total pressure 
Figure 4.6- Concluded. 
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Figure 5.1- Free vortex sheet formulation. 
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Figure 5.2- Spanwise pressure distributions. A = 63°, 

c t /c p = 0.4, a = 20°, = 0.3. From Luckring, et al. (1982). 
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Figure 5.3- Calculation of vortex core trajectory. From 
Luckring, et al. (1982). 
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Figure 5.4- Calculation of force/moment properties. 
AR = 1 delta wing. 
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a) vortex geometry 

Figure 6.1- Sheet rollup effects on global solution properties. 
From Johnson, et al . (1980). 
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Figure 6.1- Concluded. 
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Figure 6.2- Comparison of vortex sheet trajectories. 
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a) angular velocity 

Figure 6.3- Inner- to outer-flow matching 
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b) axial velocity 
Figure 6.3- Continued 
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c) circulation gradient 
Figure 6.3- Concluded. 



Figure 6.4- Nonconical effects. < 1. 
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Figure 6.5- Representative boundary conditions. A - 75°, 
a = 15°. 
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Figure 6.6- Longitudinal scaling effect. 
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a) radial velocity 

Figure 7.1- Representative delta wing solution - scaled radial 
distribution. A = 75°, a = 15°. 
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b) circulation parameter 
Figure 7.1- Continued 
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c) angular velocity 
Figure 7.1- Continued 
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d) axial velocity 
Figure 7.1- Continued 
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e) static pressure 
Figure 7.1- Continued 
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f) total pressure 
Figure 7.1- Concluded 
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a) radial velocity 

Figure 7.2- Representative delta wing solution 


distribution. A = 75 l 


= 15°. 


unsealed radial 



b) circulation parameter 
Figure 7.2- Continued 
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c) angular velocity 
Figure 7.2- Continued 
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d) axial velocity 
Figure 7.2- Continued 




































b) pressures 
Figure 7.6- Concluded. 
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